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ABSTRACT 



Aims. With our time-dependent model atmosphere code PHOENIX, our goal is to simulate light curves and spectra of hydrodynamical 
models of all types of supemovae. In this work, we simulate near-infrared light curves of SNe la and confirm the cause of the secondary 
maximum. 

Methods. We apply a simple energy solver to compute the evolution of an SN la envelope during the free expansion phase. Included 
in the solver are energy changes due to expansion, the energy deposition of y-rays and interaction of radiation with the material. 
Results. We computed theoretical light curves of several SN la hydrodynamical models in the I, J, H, and K bands and compared 
them to the observed SN la light curves of SN 1999ee and SN 2002bo. By changing a line scattering parameter in time, we obtained 
quite reasonable fits to the observed near-infrared light curves. This is a strong hint that detailed NLTE effects in IR lines have to be 
modeled, which will be a future focus of our work. 

Conclusions. We found that IR line scattering is very important for the near-infrared SN la light curve modeling. In addition, the 
recombination of Fe III to Fe II and of Co III to Co II is responsible for the secondary maximum in the near-infrared bands. For future 
work the consideration of NLTE for all lines (including the IR subordinate lines) will be crucial. 

Key words, stars: supernovae: general - radiative transfer - methods: numerical - stars: atmospheres 



1. Introduction 

We focus on the near-infrared light curves of type la supernovae 
(SNe la). The light curves in the near-infrared are of particular 
interest, because they have been claimed to be n ear standard can- 
dles a t the time of B-band maximum in the IR dKrisciunas et al.l 
2004a). In addition, many, but not all SNe la exhibit a secondary 
maximu m in the IR bands . The secondary maxi mum was first 
noted by Elia s et all (I1981I) . iHoflich et aTT(ll995l) explained the 
secondary maximum as due to the ex pansion of the IR pseudo- 
photosphere, whereas TSuntzeffl Jl996) suggest ed that it was due 
to a g lobal shift of radiation from blue to red. iPinto & Eastman! 
(2000) suggested that the transition from Fe III to Fe II as the 
dominant ionization stage was important although they implied 
that the secondary maximum occurs when the photo sphere has 
receded into the non-radioactive center. iKasenl (12006 ) performed 
a detailed study of SN la light curves in the near-infrared and 
obtained reasonable fits to observations. He also finds that the 
secondary maximum is an effect of the ionization evolution of 
iron group elements in the expanding envelope. We apply our 
time-dependent model atmosphere code PHOENIX and investi- 
gate SNe la by comparing model light curves in the near-infrared 
to observed light curves. We will also focus on the origin of the 
secondary maximum. 

We present model light curves of SNe la in the near-infrared. 
In Sect. [2] we present the methods we used in some detail. In 
Sect. [3] we present light curves of SNe la for different hydro- 
dynamical models in the I, J, H, and K bands. We consider a 
parametrized IR line scattering in the solution of the radiative 



transfer to improve our fits to observed light curves. An inves- 
tigation of the secondary maximum that has been observed in 
near-infrared light curves is discussed in a final section. 



2. Methods 

We use our time dependent model atmosphere code PHOENIX, 
version 16, to compute model light curves of type la supernovae. 
We use th e time dependen t exten sion and the method as de- 
scribed in iJack et al.1 d2009l l201ll) . We calculate SN la model 
light curves for different hydrodynamic (explosion) models. In 
our previous work, we presented model light curves in the opti- 
cal bands with reasonable fits to observations. In this work, we 
focus on the near-infrared wavelength region and model the light 
curves in the I, J, H, and K bands. 

To compute SN la light curves, our approach is a simple 
energy solver that solves for an energy change of the material 
during a time step of the envelope evolution. The energy can 
change due to the free adiabatic expansion, the deposition of en- 
ergy from y-rays from the radioactive decay of 56 Ni and 56 Co, 
and the interaction of radiation with the material. In our calcula- 
tions, we assume homologous expansion. The y-ray deposition 
is solved with the assumption of a gray atmosphere. Changes of 
ionization state are also included in our solver. 

To determine the energy change of the material due to the in- 
teraction with the radiation, we have to solve the radiative trans- 
fer. For each time step, we solve the non-gray spherical sym- 
metric radiative transfer equation for expanding atmospheres in- 
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Fig. 1. Synthetic spectrum of a SN la at day 20 plotted with the Fig. 2. Model and observed light curves in the I band. Two dif- 
filter functions of the I, J, H, and K band. ferent explosion models were used to compute the model light 

curves. 



eluding special relativistic effects. We typically use about 2,000 
wavelength points to keep the computation time reasonable. 

Given the problems with the simple LTE treatment of line 
scattering as shown in iJack et ail d201 lb . we investigate the ef- 
fects of IR line scattering on SN la light curve s in the near- 
infrared. In our previous paper (IJack et alJl20TTh . we indicated 
how important UV/optical line scattering is during the later 
phases of the light curve evolution. Here, we study this effect 
in more detail. 

The source function of the radiative transfer equation includ- 
ing scattering for an equivalent two level atom can be written as 



ight curves in J band 



(1 - e A )J A + e A B A 



(1) 



The source function is represented by S A , the Planck function by 
B A and the mean intensity by J A . All these quantities depend on 
the wavelength A. Here, we used e A as the line thermal coupling 
parameter. For e A — 1 there is only true absorption and no line 
scattering takes place. The source function is then given by S A = 
B A (pure LTE). In reality, the thermal coupling parameter e A will 
vary over the whole wavelength range. For the calculation with 
PHOENIX, it is possible to set a wavelength independent factor 
e = e A = constant to approximate LTE line scattering over the 
whole wavelength range. 

For our approach to model light curves of SNe la, we use 
the results of hydrodynamical explosion simulations, as our ini- 
tial atmosphere structures which we then evolve in time. In this 
work, we u se the explosio n model results of the deflagration 
model W7 (lNomotolll98l . We also used one delayed detona - 
tion model, the model DD25 calcu lated bv lHoflich et all ((2002). 
Model DDI 6 ( Hofl ich et al.l 120021) which we studied in previ- 
ous work is not considered here, because its low yield of 56 Ni 
makes it too dim to account for normal SNe la. We used our 
hydrodynamic solver to obtai n detailed spect ra at certain days in 
the light curve evolution, (see lJack et alJ201ll for details). Light 
curves are calculated by convolving the appropriate filter func- 
tions with the theoretical spectra. We compare our synthetic light 
curves to observed light curves of two SN la events: SN 1999ee 
and SN 2002bo. The filter functions we used here are shown 
in Fig. Q] The photometric o bservations of SN 1999e e in the J 
and H band are presented in IRrisciunas et a l. (2004b). For SN 
2002bo photome tric observations in the J , H, and K band have 
been obtained by IRrisciunas et all d2004ch . 
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Fig. 3. Light curves in the J band. Two different explosion mod- 
els were used to compute the model light curves and to compare 
them to observations. 



3. Near-infrared light curves of SNe la 

We first present the results of an approach of a constant line scat- 
tering factor e = 0.8 to the light curve modeling, where we com- 
pare the results of two different hydrodynamical models. Figure 
[2] shows the calculated model light curves of two hydrodynam- 
ical models in the I band. They are compared to two observed 
light curves of SN 1999ee and SN 2002bo. Both model light 
curves produce the steep rise during the first days after the ex- 
plosion reasonably well. However, the model light curve rises 
further, although the observed light curves show a decline after 
maximum at around 13 days after explosion. During the later 
phase around 30 days after the explosion, the model light curves 
are much too bright compared to the observed light curves. There 
is also no indication of a secondary maximum in the model 
light curves. In Fig. [3] the model light curves in the J band are 
shown together with the observed light curves. Again, the first 
phase is well represented by the model light curves. The model 
light curves rise after a first maximum to a secondary maximum, 
which is much brighter and earlier than the secondary maximum 
in the observed light curves. Furthermore, the secondary max- 
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Fig. 4. Light curves in the H band. Two different explosion mod- Fig. 6. Light curves in the I band computed with different values 



els were used to compute the model light curves. 



of the line scattering factor e. 
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Fig. 5. Light curves in the K band. The model light curves of W7 
and DD25 are significantly brighter than the observations during 
the later phase. 



imum of the DD25 model is earlier than the one of hydrody- 
namical model W7. However, the difference of 2.5 mag around 
day 35 between models and observations is enormous. A similar 
picture is seen in the H and K band as shown in Fig. [4] and Fig. 


So far, we used for the calculations of our LTE model light 
curves a constant time independent factor of e = 0.8. We now use 
this LTE line scattering factor to approximate the effect of line 
scattering on the shape of the light curves in the near-infrared. 
We calculated several model light curves each with a different 
LTE line scattering factor, e. At each point of the light curve, we 
let the atmosphere structure adapt to the new conditions with its e 
until a radiative equilibrium state is reached. This led to different 
resulting atmosphere structures or more precisely temperature 
structures at each day and thus to differences in the shapes of the 
model light curves. 

Figure [6] shows the I band light curves obtained with dif- 
ferent values of e. With a smaller e, the light curve becomes 
fainter. A look at the temperature structure in Fig.|7]reveals that 
decreased e leads to a hotter atmosphere. In the hotter atmo- 
sphere, the iron and cobalt are mostly doubly ionized. As the 



Fig. 7. Temperature structure of the envelope at day 30. A 
smaller line scattering factor e leads to a hotter atmosphere. 



atmosphere cools down, the iron and cobalt recombine to Fe II 
and Co II, and the flux increases due to numerous lines of these 
species in the near-infrared. This recombinatio n effect in the 
nea r-infrared light curve has also been found bv lKa^(l2006h 
and iHoflich et alJ (ll995l). The intersection of the model light 
curves with the observed light curves shows that the best value 
of e changes during the evolution of the light curve, e decreases 
with time (i.e., line scattering becomes more and more impor- 
tant) in order to fit the observed light curves. This is reasonable, 
because the ongoing expansion leads to lower densities at later 
times, and line scattering becomes more important because of 
the thinner atmosphere. However, during the earl y phase near 
maxim um light, e has to be closer to pure LTE dBaron et al.l 
Il996h . 

We now seek to fit the entire light curve with a single func- 
tion for e, that is we choose e = eof(t) where /(f) is a function of 
time. Our model light curves now reproduce the observed light 
curves much more faithfully. Even the secondary maximum is 
well reproduced for each band. For instance, we show in Figure 
[12] the new best fitting light curve in the I band. The solid line 
shows the W7 model light curve, which has a an e that decreases 
with time. We chose the light curve of SN 1999ee to obtain the 
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Fig. 8. Spectral evolution in the I band from day 20 to day 50 Fig. 9. Spectral evolution in the J band from day 20 to day 45 
after the explosion. The solid line displays the filter function. after the explosion. The solid line shows the filter function. 



best fit. It is clear that different values of /(f) can be used to 
obtain a better fit to the light curve of SN 2002bo. 



4. Secondary maximum in the near-infrared SN la 
light curves 

In all bands in the near-infrared, the observed light curves show a 
secondary maximum around 30 to 40 days after explosion. Since 
with a time-dependent e our model light curves reproduce the 
observed light curves pretty well, we will focus on the specific 
cause of the secondary maximum in each IR band. We investi- 
gate the cause of the secondary maximum in each of the near- 
infrared bands by looking at the spectral evolution of the model 
light curv es in the resp ective band. Basically, we confirm the 
results of lKasenl |2006). However, we can additionally assign 
features of certain elements to the individual bands and show 
detailed spectra. 

The I band is known to be dominated by the Ca IR triplet. In 
Fig. [8] the spectral evolution in the I band is shown. The broad 
feature at 8500 A is the Ca triplet. This feature stays present as 
the flux outside of the feature declines while the SN la envelope 
evolves. This leads to the decline in the light curve after the first 
maximum. Starting at day 40 a new "W-shaped" feature emerges 
at around 7500 A. These are lines of Fe II, and the feature is even 
more significant in the spectrum at day 50. From day 30 to day 
50 the temperature in the expanding envelope is decreasing. The 
temperature eventually reaches the temperature/pressure regime 
where Fe III recombines to Fe II. This increases the brightness 
in the I band due to emission lines of Fe II present in this wave- 
length region. Thus, emission lines of Fe II are responsible for 
the rise to a secondary maximum in the I band. 

We now examine the secondary maximum in the J band. The 
spectral evolution in that wavelength range is shown in Fig. [9] 
With the SN la evolving, the flux is decreasing in the J band 
because the temperature of the envelope decreases. At day 30, 
new features arise, which increase the brightness in the J band. 
These features are caused by lines of Fe II and Co II. Therefore, 
the recombination of Fe III to Fe II is responsible for the rise in 
the J band. Additionally, Co III recombines to Co II at about the 
same temperature and therefore time in the light curve evolution. 
Thus, Co II lines emerge in the J band due to this recombination. 
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Fig. 10. Spectral evolution in the H band from day 16 to day 45 
after the explosion. The solid line shows the filter function. 



The secondary maximum in the J band is caused by both emerg- 
ing Fe II and Co II lines. 

The spectral evolution in the H band is shown in Fig.[T0] In 
the spectrum at day 26, two broad features emerge in this wave- 
length range. In the spectrum of day 35 and day 45, these two 
features are clearly visible. They cause the rise in the brightness 
of the H band to a secondary maximum. We found that these 
features are caused by lines of Co II. Thus, the recombination of 
Co III to Co II causes the rise to a secondary maximum in the 
light curve of the H band. 

In the K band, the spectral evolution of our SN la model is 
shown in Fig.QT] The flux is decreasing as the envelope expands 
and cools down. However, at day 40 and day 45 features emerge 
and cause a rise in the brightness of the K band. These features 
are again lines of Co II. The recombination of Co III to Co II 
causes the secondary maximum in the K band. 

4. 1 . Non-radioactive core 

The W 7 model has a non-radioactive iron core. lPinto~& Eastman 
(2000) showed that the consequent negative temperature gradi- 
ent leads to a further reduction in the mean opacity. This may 
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Fig. 11. Spectral evolution in the K band from day 16 to day 45 
after the explosion. The solid line shows the filter function. 
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Fig. 12. Light curve in the I band with and without the non- 
radioactive iron core. 



have an influence on the duration or magnitude of the IR sec- 
ondary maximum. Therefore, it is worthwhile checking wether 
the non-radioactive center affects the shape of the secondary 
maximum in the near-infrared. We replaced the iron core of the 
W7 model with radioactive 56 Ni and computed the light curve. 
This increases the total mass of 56 Ni from 0.56 to 0.66 solar 
masses. As shown in Fig. [12] there are only small difference be- 
tween the light curves with and without the non-radioactive iron 
core. This simple variation does not capture the full effect of 
the presence or absence of a non-radioactive core, since a shift 
in the nickel distribution can have important effects on the light 
curve (E. Baron et al., in preparation). Additionally, late time 
spectra show that tha t a non-radioactive core seems to be re- 
quired in most SNe la dHoflich et al.l2004lMotohara et aLl2006t 
Fesenetal ] l2007tlGerardv et alJl2007t iMaeda et alJ2010h . 

5. Conclusion 

We applied our time-dependent model atmosphere code 
PHOENIX to model light curves of type la supernovae in the near- 
infrared wavelength range. In a first approach with a constant 
LTE line scattering parameter e in the radiative transfer, we re- 
produced the observed light curves during the first phase quite 



well. However, the model light curves during the later phase 
were too bright, and we could not reproduce the secondary max- 
imum. 

We show that a more detailed treatment of IR line scattering 
is very important for the modeling of the later phase of the near- 
infrared light curves of SNe la. We used as an approximation 
an LTE line scattering parameter that decreases in time. This 
is a good approach, because the atmosphere becomes thinner 
as expansion goes on and scattering becomes more important. 
We use this approximation to obtain fits to the observed light 
curves of SN 1999ee and SN 2002bo, improving our model light 
curves significantly. The next step is to treat the atmosphere in 
full non-LTE, where the temperature needs to adapt to the non- 
LTE condition, which requ ires substantial time on parallel su- 
percomputers. |Kai3(2006) also found that his models overshot 
the observed secondary maximum in H and K, while he mainly 
attributed this to incomplete line lists in the IR, he also suggested 
that his treatment of LTE might be a possible cause. 

The secondary maximum in each band of the near-infrared 
was quite accurately reproduced by our model light curves with 
time-dependent e. We investigated the spectral evolution and 
found that the I band secondary maximum arises due to recom- 
bination of Fe III to Fe II. An Fe II feature emerges in this band 
and increases the brightness. For the secondary maximum in the 
J band, a mix of Fe II and Co II lines emerges and causes the rise 
to a secondary maximum in the model light curve. The bands 
H and K also show a secondary maximum. This is caused by 
the recombination of Co III to Co II. We confirmed that ion- 
ization stage changes in iron group elements a re responsible for 
the secondary maximum as previously found (Kasen 2006|). In 
fact, since we use neither expansion opacities, nor the Sobolev 
approximation, our results robustly indicate the need for a full 
NLTE treatment. Furthermore, we can explicitly assign lines of 
different elements to different bands. These results also show 
how important the treatment of Fe and Co in non-LTE is and 
should be used in future work, when faster computers are avail- 
able. 
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